High-Frequency Impedance of Driven Superlattices 
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The complex impedance of a semiconductor superlattice biased into the regime of negative 
differential conductivity and driven by an additional GHz ac voltage is computed. From a simulation 
of the nonlinear spatio-temporal dynamics of traveling field domains we obtain strong variations 
of the amplitude and phase of the impedance with increasing driving frequency. These serve as 
fingerprints of the underlying quasiperiodic or frequency locking behavior. An anomalous phase 
shift appears as a result of phase synchronization of the traveling domains. If the imaginary part of 
the impedance is compensated by an external inductor, both the frequency and the intensity of the 
oscillations strongly increase. 
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Semiconductor superlattices (SL) show pronounced 
negative differential conductivity (NDC)EJ. If the total 
bias is chosen such that the average electric field is in 
the NDC region, stable inhomogeneous field distribu- 
tions (field domains)B or self-sustained oscillationscLwith 
frequencies up to 150 GHz at room temperatureQ ap- 
pear. Which of these scenarios occurs depends on bias, 
doping—temperature, and the properties o£-the injecting 
contactOiLl For a recent overview see Ref.lill. 

In order to apply the self-sustained oscillations as a 
high-frequency generator in an electronic device, it is 
crucial to know the response of the SL in an external 
circuit. A key ingredient for the analysis is the com- 
plex impedance of the SL in the respective frequency 
range. This is the subject of this paper where the com- 
plex impedance is evaluated numerically by imposing an 
additional ac bias to the SL. The interplay of the self- 
sustained oscillations and the external frequency causes 
a variety of interesting phenomena such as frequency 
locking, quasi-periodic and chaotic behavior ptbjch has 
been extensively—studied both theoreticallycj~tj and 
experimentallyiL3't3. In contrast to those studies we fo- 
cus on the response to the circuit and concurrent phase 
synchronization phenomena in this work. 

We describe the dynamical evolution of the SL by 
rate equations for the electron densities in the quan- 
tum wells, together with Poisson's equation for the elec- 
tric fields. The current densities Jj^j+i between adja- 
cent quantum wells are evaluated within the model of 
sequential tunneling, and Ohmic boundary conditions 
with a contact conductivity a are used. For details 
see Refs.EjO. Here we apply a periodic bias signal 
U(t) = Udc + U ac sin(27r^ii) and study the total current 
I(t) = AJ2j= Jj^j+i/( N + !)■ In particular we con- 
sider the Fourier component 7 ac (i^i) sin(27Wit — <fr) which 
gives the complex impedance 

As an example we consider the SL structure studied in 



Ref.Q. It consists of N = 120 GaAs wells of width 4.9 
nm separated by 1.3 nm AlAs barriers. The sample is 
n-doped with a density of 8.7 x 10 10 cm~ 2 per period and 
the sample cross section is A — 64(^m) 2 . We estimate 
the energy broadening r — 15 meV, which effectively 
gives the sum of phonon, impurity and interface rough- 
ness scattering rates. All calculations are performed at 
room temperature. 

In Fig. [l]a we display the calculated current-voltage 
characteristic for C/ ac = 0. While for large values of the 
contact conductivity a stationary field domains occur, 
self-sustained current oscillations (where pairs of accu- 
mulation and depletion fronts travel through the sam- 
ple) are found for lower a (Fig. |l|b). This_acenario is 
quite general for moderately doped sampleslZEd. In the 
following we use a — 20 A/Vm giving best agreement 
with the measured current-voltage characteristic. 

At £/dc = 2 V, the self-sustained oscillations exhibit a 
frequency of vq — 2.1 GHz. Now we study the change in 
the current signal by imposing an additional ac-bias am- 
plitude Unc with frequency v\. Fig. displays the abso- 
lute value and the phase of the complex impedance given 
by Eq. ([!]). For both strengths of the ac amplitude one 
observes strong variations in the amplitude and phase of 
Z. In particular one observes pronounced minima in \Z\ 
and a corresponding monotonic increase in (f> around fre- 
quencies v\ which are approximately-, integer multiples of 
vq . This is due to frequency lockingll3, where the external 
frequency v\ modifies the main oscillation frequency to 
Vq such that v\/vq is a rational number. While locking 
occurs in a rather wide frequency range for v\/v^ 6 N, 
locking into rational numbers with larger denominators 
are less easy to observe (The corresponding widths of 
the locking intervals roughly follow th£_.position in the 
Farey-tree of the rational numbers, seeEZl. The width of 
the locking ranges also increases with t7 ac .) A more de- 
tailed examination indicates that the local minima in \Z\ 
correspond to the locking regions, and <f> increases with 
V\ in those intervals, while a decrease is frequently ob- 
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served outside these ranges where quasi-periodic or possi- 
bly chaotic beha*ipr occurs. Thus the presence of several 
locking intervalst3 explains the variety of peaks. 

For high frequencies v\ ^> v$ the variations of Z be- 
come less pronounced and \Z\ ~ 300S1 and <f> « — 7r/2 is 
observed. This corresponds to a capacitive current with 
an effective frequency-dependent capacitance C e ff(fi) ~ 
1/{2-kvi\Z\). At v x = 10 GHz this gives C cS « 50 fF. 
This capacitance results from the interaction with the 
internal front dynamics in the SL structure; the intrinsic 
sample capacitance e r eoA/ (Nd) = 10 fF is significantly 
smaller. 

While we have only shown results for Udc — 2 V and 
a = 20 A/Vm here, we checked that the features dis- 
cussed above neither change for different biases nor dif- 
ferent contact conductivities, albeit the main oscillation 
frequency Vq and the locking ranges change slightly. 

Now we focus on the 1/1 locking region which is of 
interest if the device is used as an oscillator. In Fig. || 
we show the current and bias signal for different frequen- 
cies v\. In parts (a) (Vi = 1.95 GHz) and (f) {v\ = 2.3 
GHz) no complete locking occurs and the current does 
not exhibit a periodic signal. Between these frequencies 
the current signal is periodic with a frequency = v\ 
imposed by the external bias. While the current peaks 
occur around the minima of the external bias for v\ = 2 
GHz (corresponding to (j> w — 7r), the delay between the 
current peaks and the bias maxima decreases with fre- 
quency until they are approximately in phase at the up- 
per boundary of the locking range at v\ = 2.25 GHz 
(corresponding to <p w 0). It is intriguing to note that 
this behavior is opposite to the response of a damped lin- 
ear oscillator, where the phase between the driving sig- 
nal and the response shifts from to n with increasing 
driving frequency. The current peaks correspond to the 
formation of a domain at the emitter, while during the 
current minima the domain traverses the sample. The 
domain transit velocity is larger for smaller voltage. The 
phase shift between current and voltage adjusts such that 
the domain velocity increases with increasing driving fre- 
quency during the whole locking interval. This is a phase 
synchronization effect induced by the domain dynamics. 

Let us analyze the behavior close to the onset of the 
locking interval in detail, see Fig. § The main frequency 
component of the current signal is given by z/q. In the 
locking region it is equal to the frequency of the driving 
bias v\ such that i/q — v\ vanishes in a finite range of v\ . In 
contrast, far away from the locking is essentially given 
by the free oscillation frequency vo, thus Uq — v\ exhibits 
a linear relation. Close to the boundaries o f the lock- 
ing region a square root behavior v$ — v\ cx \J\v\— z/J rit | 
can be detected, p-which is a general feature of frequency 
locking see, e.g.Jla. 

Fig. [|b shows the amplitude and the monotonically 
increasing phase of the impedance over the whole lock- 
ing interval. We extract an impedance Z — il50il for 
vo = 2.1GHz. Now we compensate this impedance by 
an external circuit with an inductor of L = 10 nH and a 



resistor of 25 f2 in series with the SL. (The resistor may 
result both from connecting wires or due to radiation 
damping in a real device.) Fig. |5| shows that the oscilla- 
tion mode is strongly affected: The frequency is increased 
to v = 5.6 GHz, the oscillation is more sinusoidal, and 
the current amplitude increases. Thus the device per- 
formance of the SL is strongly improved. This effect is 
due to a different oscillation mode, where the domains 
are quenched such that they only traverse a small part 
of the SL. As the domain velocity is almost constant, the 
frequency increases. 

We have shown that the amplitude and the phase of 
the complex impedance Z exhibit strong variations with 
the driving frequency v\ . The distinct minima of the am- 
plitude of the impedance correspond directly to the fre- 
quency locking intervals. In these regions the frequency 
components of the current at v\ are particularly large, 
leading to small \Z\. The phase of the impedance is 
monotonically increasing with frequency from — 7r to zero 
in the locking intervals, which constitutes a phase shift 
of 7r compared to the standard behavior of a driven oscil- 
lator. These effects are caused by phase synchronization 
of the spatio-temporal dynamics of the traveling field do- 
mains in the NDC regime. Compensating the imaginary 
part of the impedance in an external circuit strongly im- 
proves the device performance of the SL. 
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FIG. 1. (a) Calculated current-voltage characteristics for 
different contact conductivities cr, and £/ ac = 0. For a = 20 
A/Vm and 40 A/Vm self-sustained current oscillations are 
found in the indicated bias ranges. In this case we display 
the time-averaged current, (b) Self-sustained oscillations for 
(7dc — 2 V, (j = 20 A/Vm. The time series of the current I(t) 
and the space-time plot of the electron densities ni(t) in the 
quantum wells is shown. Black indicates low electron density 
(depletion front), white indicates high electron density (accu- 
mulation front). The emitter is at the bottom, the collector 
at the top. 



IUUU 



(a) V, = 1 .95 GHz 



(c) V, = 2.05 GHz 




(e) V, = 2.25 GHz 




20 10 
15 8 



5> 5? I 
10" E 




10 — E 



(d) v,= 2.15 GHz 




(f) V1 = 2.30 GHz 
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FIG. 3. Current / (full line) and driving bias U 
(dashed line) versus time for different driving frequencies v\ 
(L/ac = 0.4V, U dc = 2V). 
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FIG. 4. (a) Shift of the fundamental frequency Vq of the 
current signal and (b) amplitude and phase of the complex 
impedance as a function of the driving frequency v\ in the 
vicinity of the 1/1 locking region (J7 a c = 0.2 V, Udc = 2V). 




v 1 [GHz] 



FIG. 2. Amplitude \Z\ (full line) and phase (f> (dashed line) 
of the complex impedance Z as a function of the driving fre- 
quency for two different values of the ac bias: (a) U ac = 0.2V, 
(b) U ac — 0.6V. The shaded areas indicate locking intervals 
marked by v^jvl (U dc = 2V, a = 20 A/Vm). 
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FIG. 5. Operation of the SL device with an inductor of 10 
nH and a resistor of 25 Q in series (Udc = 2V, (V ac = 0) (a) 
Current versus time and space-time plot of the electron den- 
sity in the SL (as in Fig. lb), (b) Electric field versus position 
for different times (as indicated in a). 
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